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Results of neutrino oscillation experiments have always been presented on the (sin 2 28, Am 2 ) param- 
eter space for the case of two-flavor oscillations. We point out, however, that this parameterization 
misses the half of the parameter space J < < § ( "the dark side" ) , which is physically inequivalent 
to the region < 6 < j ("the light side") in the presence of matter effects. The MSW solutions 
to the solar neutrino problem can extend to the dark side, especially if we take the conservative 
attitude to allow higher confidence levels, ignore some of the experimental results in the fits, or relax 
theoretical predictions. Furthermore, even the so-called "vacuum oscillation" solution distinguishes 
the dark and the light sides. We urge experimental collaborations to present their results on the 
entire parameter space. 



In the Standard Model of particle physics, neutrinos 
are strictly massless. Recently, however, the Super- 
Kamiokande collaboration studied atmospheric neutrinos 
and reported a strong evidence for neutrino oscillations 
0, and hence a finite neutrino mass. The most likely 
interpretation of their data is the oscillation between 
and v T . This made it also natural to interpret another 
long-standing issue in neutrino physics, the deficit of the 
solar v e flux in terms of neutrino oscillations. How- 
ever, the solar neutrino deficit has not been regarded as 
convincing evidence for neutrino oscillations in the com- 
munity. The reason is probably multifold but two main 
objections are the following. (1) Neutrino experiments 
are so difficult that it is possible that some of the data 
are not entirely correct. (2) The physics of the Sun is so 
complex that the neutrino flux calculations in the Stan- 
dard Solar Model (SSM) may have underestimated the 
theoretical uncertainties. 

To resolve this situation, a new generation of solar 
neutrino experiments, such as Super-Kamiokande, SNO, 
Borexino, GNO, KamLAND, etc, is looking for an evi- 
dence for solar neutrino oscillations without relying on 
the SSM in well-understood experimental environments. 
They aim not only at establishing oscillations but also at 
overdetermining the solution in the next few years. Such 
data will eventually supersede data from the past exper- 
iments. It is therefore important to analyze the future 
data without too much prejudice based on the past data. 

In this letter, we point out that the study of neutrino 
oscillations on the (Am 2 , sin 2 20) parameter space done 
traditionally is incomplete, since it covers only the range 
< 9 < j ("the light side"). Indeed, some of the solu- 
tions to the solar neutrino puzzle extend to the other half 
of the parameter space j < 9 < ~, which we call "the 
dark side," and hence it is phenomenologically necessary 
to include both halves of the parameter space. This is 
especially true once one employs a more conservative atti- 



tude which either allows higher confidence levels, ignores 
some of the experimental data (especially Homestake j| ) , 
or relaxes the theoretical prediction on the 8 B flux. 

Neutrino oscillations occur if neutrino mass eigenstates 
are different from neutrino weak eigenstates. Assuming 
that only two neutrino states mix, the relation between 
mass eigenstates (y\ and v^) and flavor eigenstates (for 
example v e and v^) is simply given by 

^i) = cos0|^ e ) — sin#|f M ), 

\h> 2 ) = sin% e ) + cos0|i/ M ), (1) 

where 9 is the vacuum mixing angle. The mass-squared 
difference is defined as Am 2 = m| — m 2 . We are in- 
terested in the range of parameters that encompasses all 
physically different situations. First, observe that Eq. (jl]) 
is invariant under 9 — > 9 + n, \v e ) — > — |i/ e ), \v^) — ► — \v^), 
and hence the ranges [—§,§] and [-|, ^L] are physically 
equivalent. Next, note that it is also invariant under 
9 — > — 9, \ffj) — > — \v%) — > — |f2), hence it is sufficient 
to only consider 9 G [0, J]. Finally, it can also be made 
invariant under 9 — » ~ — 9, \v^) —> by relabeling 

the mass eigenstates \vi) <-» ^2), *- e - Am 2 — ► —Am 2 . 
Thus, we can take (Am 2 > 0) without loss of general- 
ity. All physically different situations are obtained by 
allowing < 6 < § . 

For the case of oscillations in the vacuum, the survival 
probability is given by 

P(y e ->• v e ) = 1 - sin 2 20 sin 2 (l.27^-l}j . (2) 

Here, Am 2 is given in eV 2 /c 4 , E in GeV, and L in km. 
In this case the oscillation phenomenon can be param- 
eterized by Am 2 and sin 2 29, since 9 and ? — 9 yield 
identical physics. Therefore we can restrict ourselves to 
< 9 < j, and use the parameter space (Am 2 , sin 2 29) 
without any ambiguity. This is indeed an adequate pa- 
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rameterization for reactor antineutrino oscillation exper- 
iments, short-baseline accelerator neutrino oscillation ex- 
periments, and v ^ «-> v T atmospheric neutrino oscillation 
experiments. 

In certain cases, however, neutrino-matter interactions 
can dramatically change the oscillation probability |Q. 
These matter effects are particularly important in ex- 
plaining the solar v e flux deficit in terms of neutrino os- 
cillations. In the presence of matter effects, Eq. (0) is 
modified to 

P(y e -> v e ) = Pi cos 2 6+ (I- Pi) sin 2 

, / A777 2 \ 

— v/P c (l - Pcj cos 20 M sin 20 cos f 2.54—— -L + 5J , (3) 

where P c is the hopping probability, 0^ is the mix- 
ing angle at the production point, Pi = P c sin 2 9 m + 
(1 — P c ) cos 2 0m, and S is a phase induced by the mat- 
ter effects, which is not important for our purposes. 
See Ref. Jj^H for notation. Because Pi depends on 
Am 2 cos 20, the half of the parameter space tradition- 
ally considered < 9 < j (the light side) is physically 
inequivalent to the other half f < < § (the dark 
side). However, all data analysis have been reported 
on the (Am 2 , sin 2 29) plane with positive Ato 2 only for 
solar neutrino experiments and hence only half of the 
parameter space has been analyzed. Even though the 
dark side has been studied in the context of three-flavor 
and four-flavor Q neutrino oscillations, the impor- 
tance of studying both halves for the simplest case of 
two-flavor oscillations has been largely ignored in the lit- 
erature. There also appeared to be a misconception in 
the literature that physics was discontinuous at maximal 
mixing 9 = j. For instance, matter effects in the Earth 
were once thought to disappear as the mixing approached 
maximal. However, the authors of Ref. |J emphasized 
that the matter effects remain important even for the 
maximal mixing, and the present authors further showed 
that physics is completely continuous beyond = ? 
One can still retain the dark side with only < 9 < ? 
if a separate parameter space with Am 2 < is added. 
This is indeed what Super-Kamiokande did in the case of 
v s oscillations of atmospheric neutrinos [ [To[ . How- 
ever, as argued in ||], it is more natural to use < 9 < ^ 
with the fixed sign of Am 2 to exhibit the continuity of the 
physics between the two halves of the parameter space. 

Part of the reason why the dark side has been neglected 
in the literature is that it is impossible to obtain v e sur- 
vival probabilities less than one half when the two mass 
eigenstates are incoherent, i.e., when the last term in 
Eq. (|) is absent. (This occurs in the so-called "MSW 
region" 1CT 8 < Am 2 < 1CT 3 eV 2 pj.) Indeed, the data 
from the Homestake experiment p| used to be about a 
quarter of the SSM prediction, and this could have been 
an argument for dropping the dark side entirely in the 
MSW region. However, the change from BP95 Jl^] to 



BP98 (T3) calculations increased the Homestake result to 
about a third of the SSM with a relatively large theoret- 
ical uncertainty. Therefore it is quite possible that the 
"MSW solutions" extend to the dark side as well. More- 
over, some people question the SSM and/or the Homes- 
take experiment, and perform fits by ignoring either (or 
both) of them G3j. We show below that some of the 
MSW solutions indeed extend to the dark side and hence 
it is necessary to explore the dark side experimentally. 
If we further relax the theoretical prediction on the 8 B 
solar neutrino flux and/or ignore one of the solar neu- 
trino experiments in the global fit, the preferred regions 
extend even deeper into the dark side. 

Another possible reason for disregarding the dark side 
is that the so-called "vacuum oscillation region" (Am 2 < 
10~ 9 eV 2 ) was believed to be the same in the light 
and dark sides. This is because Pi approaches cos 2 9 
for Am 2 < 1CT 9 eV 2 (P/MeV) and Eq. §) reduces 
to Eq. (j|). It is remarkable, however, that low-energy 
(especially pp) neutrinos do not reach this limit for 
Am 2 > 10~ 10 cV 2 and hence the preferred regions are 
different in the light and the dark sides This ob- 

servation also implies that the separation of the MSW 
region and the vacuum oscillation region as traditionally 
done in the global fits is artificial and misleading. It is 
important to study the entire range of Ato 2 continuously. 

If sin 2 29 is not a good parameter, what is the alterna- 
tive? Two suggestions have been made in the literature. 
One is sin 2 0, which is natural since the matter effect de- 
pends directly on sin 2 || . If plotted on the linear scale, 
pure vacuum oscillations would yield physics reflection- 
symmetric around sin 2 = 0.5. If plotted on the log 
scale, the reflection symmetry is lost, but it is still a useful 
parameterization as physics is completely continuous and 
smooth from the light to the dark side. Another possible 
parameterization is tan 2 0, which retains the reflection 
symmetry for pure vacuum oscillation around tan 2 = 1 
if plotted on the log scale . We employ tan 2 for the 
analysis below because we would like to use the log scale 
to present the MSW solutions as well as the importance 
of the matter effect on the "vacuum oscillation" region 
at the same time. Note that the Jacobian from sin 2 
or tan 2 9 to sin 2 20 is singular at = j and plots with 
sin 2 20 will display unphysical singular behavior there J| . 

We next present the results of global fits to the cur- 
rent solar neutrino data from water Cherenkov detectors 
(Kamiokande and Super-Kamiokande) a chlorine 

target (Homestake) H and gallium targets (GALLEX 
and SAGE) jl7| on the full parameter space. We do not 
include the spectral data from Super-Kamiokande [|l8| as 
it appears to be still evolving with time. The fit is to 
the event rates measured at these experiments only. In 
computing the rates we include not only the pp, 7 Be, 
and 8 B neutrinos, but also the 13 N, 15 0, and pep neutri- 
nos. We use Eq. (0) with P c computed in the exponential 
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approximation for the electron number density profile 
in the Sun, and properly account for neutrino interac- 
tions in the Earth during the night with a realistic Earth 
electron number density profile by numerically solving 
Schrodinger equation as described in |l. Since the mix- 
ing angle at the production point in the Sun's core de- 
pends on the electron number density, we integrate over 
the production region numerically. We treat the corre- 
lations between the theoretical uncertainties at different 
experiments following Ref. To insure a smooth tran- 
sition between the MSW and the vacuum oscillation re- 
gion, we integrate over the energy spectrum (including 
the thermal broadening of the 7 Be neutrino "line") for 
Am 2 < 10~ 8 eV 2 and average the neutrino fluxes over 
the seasons. For Am 2 > 10 -8 eV 2 we treat the two mass 
eigenstates as incoherent. Results are completely smooth 
at Am 2 = 10~ 8 eV 2 , as expected. This allows us to fit 
the data from Am 2 = fCP n -fO~ 3 eV 2 all at once, un- 
like previous analyses which separate out the "vacuum 
oscillation region" from the rest. 

As was mentioned earlier, we take the global fit to the 
currently available data only as indicative of the ulti- 
mate result because we expect much better data to be 
collected in the near future to eventually supersede the 
current data set. We would like to keep our minds open to 
surprises such as the possibility that one of the earlier ex- 
periments was not entirely correct or that the theoretical 
uncertainty in the flux prediction was underestimated. 
In this spirit, we employ more conservative attitudes in 
the global fit than most of the analyses in the literature 
in the following three possible ways. (1) We allow higher 
confidence levels, such as 3 a. (2) We relax the theoret- 
ical prediction on the neutrino flux. (3) We ignore some 
of the experimental data in the fit. 

The global fit results are presented in Fig. at the 2 a 
(95% CL) and 3 a (99.7% CL) levels defined by x 2 ~xL™ 
for two degrees of freedom. It is noteworthy that both 
the LMA and LOW solutions (we use the nomenclature 
introduced in ]l9[ ) extend to the dark side at the 3 a 
level. At 99% CL, however, the LMA solution is confined 
to the light side. This result is consistent with the two- 
flavor limit of the three-flavor analysis in pj and the four- 
flavor analysis in || , where the spectral data is included 
and the LOW solution extends into the dark side at 99% 
CL. Another interesting fact is that the LOW solution 
is smoothly connected to the VAC solution, where the 
preferred region is clearly asymmetric between the light 
and the dark sides. Note that, at Am 2 ~ 1CP 9 eV 2 , 
the allowed region is bigger in the dark side. The region 
10~ 9 < Am 2 < 1CT 8 eV 2 was, to the best of the authors' 
knowledge, never studied fully in the literature and this 
result demonstrates the need to study the entire Am 2 
region continuously without the artificial separation of 
the "MSW region" and "vacuum oscillation region," as 
traditionally done in the literature. 

We next present a fit where the theoretical prediction 
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FIG. 1. A global fit to the solar neutrino event rates at 
chlorine, gallium and water Cherenkov experiments. The re- 
gions are shown at 2 a (light shade) and 3 a (dark shade) 
levels. The region tan 2 6 > 1 is the dark side 8 > j. 

of the 8 B flux is relaxed. Even though the helioseismol- 
ogy data constraints the sound speed down to about 5% 
of the solar radius |l3| , the core region where 8 B neutri- 
nos are produced is still not constrained directly. Given 
the sensitive dependence of the 8 B flux calculation on the 
core temperature $s B ~ T 22 , we may consider it as a free 
parameter in the fit. This can be done within the formal- 
ism of Ref. |7j by formally sending the error in Cb c to 
infinity. The result is presented in Fig. ||. The preferred 
region extends more into the dark side than the previ- 
ous fit. Even though the LMA and LOW solutions are 
connected in this plot, the lack of a large day-night asym- 
metry at Super-Kamiokande would eliminate the range 
3 x 1(T 7 < Am 2 < KT 5 eV 2 for 0.2 < tan 2 9 < 1 @. 
It is important for Super-Kamiokande to report their ex- 
clusion region on the dark side. 

Finally, we present a fit where the event rate measured 
at the Homestake experiment is not used in Fig. |[ This 
may be a sensible exercise given that the neutrino capture 
efficiency was never calibrated in this experiment. The 
preferred region extends into the dark side even at the 
95% CL. Note also the asymmetry between the dark and 
the light sides even for Am 2 < 10 -9 eV 2 . 

We expect the data of the current and next gen- 
eration of solar neutrino experiments, such as Super- 
Kamiokande, SNO, GNO, Borexino, KamLAND, to 
eventually supersede the current data set. Therefore we 
regard the above global fits only as estimates of the ulti- 
mate results. The most important point is that all exper- 
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FIG. 2. A global fit to the solar neutrino event rates at 
chlorine, gallium and water Cherenkov experiments, where 
the 8 B flux is treated as a free parameter. Contours are shown 
at 2 a (light shade) and 3 a (dark shade). 
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FIG. 3. A global fit to the solar neutrino event rates at 
the gallium and water Cherenkov experiments but not at the 
chlorine experiment. Contours are shown at 2 a (light shade) 
and 3 a (dark shade). 



imental collaborations should report their results, both 
exclusion and measurements, on both sides of the param- 
eter space, without unnecessary theoretical bias towards 
the light side. We strongly urge the experimental collab- 
orations to consider this point. 
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